Off-resonant entanglement generation in a lossy cavity 
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We provide an exact and complete characterization of the entanglement dynamics of two qubits 
coupled to a common structured reservoir at zero temperature. We derive the conditions to max- 
imize reservoir-induced entanglement for an initially factorized state of the two-qubit system. In 
particular, when the two qubits are placed inside a lossy cavity, we show that high values of entangle- 
ment can be obtained, even in the bad cavity limit, in the dispersive regime. Finally we show that, 
under certain conditions, the entanglement dynamics exhibits quantum beats and we explain their 
physical origin in terms of the interference between two different transitions coupling the dressed 
states of the system. 

PACS numbers: 03.67.Bg, 03.65.Yz 



I. INTRODUCTION 

Quantum entanglement is the powerful resource lying 
at the root of a new class of technologies based on the 
laws of quantum theory. The coherent manipulation of 
quantum systems involves very delicate procedures since 
the inevitable interaction with their surroundings leads 
to a loss of information that causes both the transforma- 
tion of quantum superpositions into statistical mixtures, 
a process called decohcrence, and the disappearance of 
quantum entanglement in composite systems. 

Recently, it has been shown that entanglement can be 
lost completely in a finite time despite the fact that com- 
plete decoherence only occurs asymptotically. This phe- 
nomenon, named entanglement sudden death, has been 
theoretically predicted by Yu and Eberly [l[ , and exper- 
imentally observed for entangled photon pairs [2j and 
atomic ensembles Q. Typically, entanglement sudden 
death occurs when the two qubits interact with two in- 
dependent environments as for the case, e.g., of two en- 
tangled qubits placed inside two different cavities. For 
such a configuration, a class of states has been identified 
which do not experience a complete entanglement loss 
despite the interaction with local vacuum environments 
However, for finite temperature environments the 
sudden death occurs almost independently of the initial 
state of the qubit pair || , although with details that can 
depend on the amount of non-Markovianity of the en- 
vironments @. In this context, a deeper understanding 
of the sudden death process has been gained by looking 
at the quantum correlations shared by the environments 
which show a sudden birth (though with a quite counter- 
intuitive timing) @, [1] . 

A completely different phenomenology emerges when 
the qubits interact with the same environment. In 
this case, indeed, entanglement can be created start- 
ing from a factorized state or it can even revive after 
a sudden death. This is due to the effective qubit- 
qubit interaction mediated by the common reservoir 
U El El El El, El EE 03- Many theoretical pa- 
pers have studied reservoir-induced entanglement in the 



Markovian regime, that is when the coupling between the 
qubits and the environment is weak enough to neglect the 
feedback of information from the reservoir into the sys- 
tem (memory less dynamics). An interesting extension 
to these approaches, that goes beyond the Born-Markov 
approximation, has been presented in Ref . [ItJ ■ 

In a recent paper, we have studied the dynamics of 
two qubits coupled to a common structured environ- 
ment using an exact approach that does not rely on 
the Born-Markov approximation [l8[. We focused on 
the case in which the qubits were identical and resonant 
with the cavity field, whose spectrum was modelled as a 
Lorentzian. In this paper we extend our analytical ap- 
proach to describe the more general situation in which 
the qubit frequencies are different and non-resonant with 
the main mode supported by the cavity. Our new an- 
alytical results allow us to characterize completely and 
exactly the entanglement dynamics for a generic initial 
two-qubit state containing one excitation. We study the 
time evolution of the entanglement and its dependence 
on several parameters, all in principle adjustable in the 
experiments: the relative coupling between the atoms 
and the cavity field, the initial amount of entanglement, 
the frequency of the qubits, the detuning from the cavity 
field, and the quality factor of the cavity. In this way we 
determine the conditions to achieve maximal reservoir- 
induced entanglement generation for an initial factorized 
state of the qubits, and to minimize the loss of entangle- 
ment for an initial entangled state. 

Depending on the matching of the qubit frequencies, 
we will distinguish two scenarios displaying different 
qualitative long time behavior. If the two qubits have 
the same transition frequency {uj\ = 0J2), a decoherence- 
free state (subradiant state) exists [ijj [2(| • Due to the 
presence of such a dark state, a non-zero asymptotic en- 
tanglement can be obtained in this case. On the other 
hand, if the two qubits have different transition frequen- 
cies (tt>i ^ 0J2), no subradiant state exists, so that the 
stationary entanglement always vanishes. For the sake 
of brevity, we refer to these two cases as subradiant and 
non-subradiant scenario, respectively. 
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One of our main results is the demonstration that high 
values of reservoir-induced entanglement can be obtained 
in the dispersive regime even in the bad cavity limit. 
In general, in this regime the dynamics of the concur- 
rence (that we employ to quantify entanglement (2lj ) is 
characterized by a quasi-regular and quasi-periodical pat- 
tern since the cavity photon is only virtually excited and 
therefore the two-qubit system is less affected by the cav- 
ity losses. Finally, in the good cavity limit, we predict 
the occurrence of quantum beats of entanglement and 
explain their physical origin. 

The paper is structured as follows. In Sec. II we 
present the microscopic Hamiltonian model, for which 
the exact analytical solution is presented in Sec. Ill, 
where we focus on the case in which the spectrum of the 
environment is Lorentzian as, e.g., for the electromag- 
netic field inside a lossy resonator. In Sec. IV and V 
we present and discuss our main results by looking at 
the entanglement dynamics in the subradiant and non- 
subradiant scenarios, respectively, for different coupling 
regimes and different initial states. Finally, Sec. VI con- 
tains summary and conclusions. 



II. THE MODEL 

We study an open quantum system consisting of two 
qubits coupled to a common zero-temperature bosonic 
reservoir in the vacuum. The Hamiltonian describing the 
total system is given by 



H — Hs + Hr + H u 



(1) 



where Hs is the Hamiltonian of the qubits system cou- 
pled, via the interaction Hamiltonian Hi nt , to the com- 
mon reservoir, whose Hamiltonian is Hr. 

The Hamiltonian for the total system, in the dipole 
and the rotating-wave approximations, can be written as 
(assuming H = 1) 



Hr 
Hint = 



(1) (1) , (2) (2) 



k 



^kb' k bk, 
(i 

ai<7+ 



(2) 
(3) 

(4) 



where b k , bk are the creation and annihilation operators 
of quanta of the reservoir, a± and uij are the inversion 
operators and transition frequency of the j-th qubit (J = 
1,2); finally Uk and ctjgk are the frequency of the mode 
k of the reservoir and its coupling strength with the j-th 
qubit. 

Here, the a's are dimensionless real coupling constants 
measuring the interaction strength of each single qubit 
with the reservoir. In particular, we assume that these 
two constants can be varied independently In the case 
of two atoms inside a cavity, e.g., this can be achieved by 



changing the relative position of the atoms in the cavity 
field standing wave. We denote with ax = {a\ + a?,) 1 / 2 
the collective coupling constant and with rj = aj/ax the 
relative interaction strength. 



A. Dynamics of the qubit system 

We assume that initially the qubit system and the 
reservoir are disentangled. We restrict ourselves to the 
case in which only one excitation is present in the system 
and the reservoir is in the vacuum. In this case the initial 
state for the whole system can be written as 



*(0)>= c 01 [l> 1 |0> 2 +c 02 |0> 1 |l) 



I R ' 



(5) 



where |0) ■ and (j = 1, 2) are the ground and excited 
states of the j-th qubit, respectively, while lOfe)^ is the 
state of the reservoir with zero excitations in the mode 
k. 

The time evolution of the total system, under the ac- 
tion of this Hamiltonian, is given by 



|*(t)) = ci(t)|l)i|0) a |0)H + C2(t)|0)i|l) 2 |0)j 



k 



c fc (t)|0)i|0> 2 |lfc) 



(6) 



where |lfe)_R is the state of the reservoir with only one 
excitation in the fc-th mode and \0)r = (£) k \0k}- 

The reduced density matrix describing the two-qubit 
systems, obtained from the density operator \^ (t)) (t)\ 
after tracing over the reservoir degrees of freedom, takes 
the form 



P(t) = 



/0 

| Cl (i)| 2 ci(t)c£(*) 

ct(t)c 2 (t) |c 2 (i)| 2 

\0 1 



C2 



I 2 / 



(7) 

The two-qubit dynamics is therefore completely charac- 
terized by the amplitudes ci j2 (i)- 

Introducing the j-qubit detuning from the mode k, 

6j? = LUj — Uk, the equations for the probability am- 
plitudes take the form 

c,-(t) = -ia^dke^Ckit), j = l,2 (8) 
k 

c k (t) = -igi[a x e- i5 ^ t c 1 {t) + a 2 e- iS ^ t c 2 {t)\ . (9) 

Formally integrating Eq. @ and inserting its solution 
into Eqs. ©, one obtains two integro-differential equa- 
tions for ci j2 (i), 



"i \9k\ e 



2 A5^(t-ti) 



Cl(*l) 



■uia 2 \gk\ e 



2^i 1) t e -4 2) *i C2 (t 1 ) 



(10) 
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c 2 (i) 



E 



oi\OL2 \9k\ e * e" 



L ci(ti) 
(11) 



In the continuum limit for the reservoir spectrum the 
sum over the modes is replaced by the integral 



dwJ(u)), 



where J(w) is the reservoir spectral density. In the fol- 
lowing wc focus on the case in which the structured reser- 
voir is the electromagnetic field inside a lossy cavity. In 
this case, the fundamental mode supported by the cavity 
displays a Lorentzian broadening due to the non-perfect 
reflectivity of the cavity mirrors. Hence the spectrum of 
the field inside the cavity can be modelled as 



w 2 



A 



7T (u-OJ c Y+\ 2 ' 



(12) 



where the weight W is proportional to the vacuum Rabi 
frequency and A is the width of the distribution and 
therefore describes the cavity losses (photon escape rate). 

We now introduce the correlation function fit — t±), 
defined as the Fourier transform of the reservoir spectral 
density J(u>), 

f(t -t x )= [ duJ^e^"-^-^, 



where uj c is the fundamental frequency of the cavity. In 
terms of the correlation function Eqs. (fT0|) - (fTTj) become 



<*(*) 



dti \a\ ci(ti) + aia 2 c 2 (ii)e 



xf{t-t 1 )e lSl{t - tl \ 



-i <5 2 i*i 



(13) 



Ca(t). The main steps for deriving the general solution 
are outlined in Appendix A. For specific forms of the 
reservoir spectral density, as the one we consider in this 
paper, it is possible to obtain simple analytic expressions 
for these coefficients. 

Before discussing the general features of the dynamics 
we notice that, when the two qubits have the same tran- 
sition frequency, lo\ = u> 2 , a subradiant, decoherence-free 
state exists, that does not decay in time. The existence 
of the subradiant state does not depend on the form of 
the spectral density and therefore on the resonance/off- 
resonance condition. Such a state takes the form 



|^_)=r a |l) 1 |0) 2 -r 1 |0) 1 H) s 



(17) 



When the two qubits have different frequencies, lo\ ^ u)2, 
there is no decoherence-free state. 

This simple consideration enables us to draw general 
conclusions about the dynamics of entanglement for long 
times. Indeed, one can observe two qualitatively different 
behaviors. In the subradiant scenario, occurring for u>i — 
u>2, a subradiant state exists and therefore that part of 
the initial entanglement stored in \ip-) will be 'trapped' 
for arbitrary long times. In the non-subradiant scenario, 
when loi ^ u) 2 , the subradiant state does not exist. Hence 
all initial entanglement will decay and is eventually lost 
for long times. 

We now derive the solution for the coefficients C\ (t) and 
c 2 (i) and study the entanglement dynamics discussing 
separately the two cases outlined above. 



B. Subradiant Scenario 

For w\ = LU2 the analytical solution for the amplitudes 
ci(t) and c 2 (i) takes a simple form, with a structure anal- 
ogous to the solution of the resonant case presented in 
Ref. Ql, 



c 2 (t) = - f dh [aiaaci^e 4521 * 1 + a 2 ! c 2 (ti)] 
Jo 

x/^-tiK^*-* 1 ', (14) 

where Sj = ujj — lo c and S21 = u>2 — wi. 

Performing the Laplace transform of Eqs. (|13[1 - iTT4)) 
yields 

sci(s) — ci(0) = - [ a\ ci(s) + ai«2 c 2 (s + i 5 2 \) ] 

x/O-ifc), (I 5 ) 

sc 2 (s) - c 2 (0) = - [aia 2 ci(s - i 5 2 i) + alc 2 (s)] 

xf(s-iS 2 ). (16) 

From the equations above one can derive the quantities 
ci(s) and c 2 (s)- Finally, inverting the Laplace transform 
one obtains a formal solution for the amplitudes c\ (t) and 



ci(t) =[r 2 2 + r\ £(t) ] Cl (0) - nr 2 [ 1 - £(t) } c 2 (0), 



(18) 



c 2 (t) = -nr 2 [ 1 - £(t) } ci(0) + [r\ + r\ £{t) } c 2 (0), 



(19) 



with 



£(t) 



-(A— iS) t/2 



where #1 



cosh(m/2) 



5 and SI 



X-iS 



sinh (fit/ 2) 



(20) 



2^2 



sj\ 2 - W R - mx, with 

S 2 the generalized Rabi frequency and 
1Z = Wax the vacuum Rabi frequency. 

As in the resonant case, the state does not evolve 
in time and the only relevant time evolution is the one of 
its orthogonal superradiant state 



|V+> = ri|l)i|0) 2 + r 2 |0)i|l) 



(21) 
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The function £{t) is the survival amplitude of the super- 
radiant state (ip + (t)\ip + (0)) = £{t). If we express the 
initial state of the qubits as a superposition of \ip±), that 
is |V(0)) = /3-|V-) +0+\1>+) with (3 ± = (V>±|V>(0)>, we 
see that, while part of the initial state will be trapped 
in the subradiant state | ?/>-), another part will decay fol- 
lowing Eq. (|20p . Thus the amount of entanglement that 
survives depends on the specific initial state and on the 
value of the coefficients Tj . 

C. Non-subradiant Scenario 

For lo\ 7^ u>2 no subradiant or decohcrence-free state 
exists and, as a consequence, the analytical expression 
for the amplitudes ci^(t) becomes more complicated 

ci(t) = £n(t;ri)d(0)+£ 12 (t;ri)ca(0), (22) 
c 2 (t) = f 21 (t;ri)ci(0)+faa(*;ri)ca(0), (23) 

where the functions £ij(t;r\) depend not only on time 
but also on the value of ri . 

We emphasize that in both scenarios, the solution of 
the differential equations for the amplitudes c\^{t) is ex- 
act as we have not performed neither the Born nor the 
Markov approximation. The structure of the functions 
£ij{t',ri) and the main steps to the solution are briefly 
outlined in Appendix A. 

D. Dispersive regime 

In this subsection we focus on the system dynamics 
when the qubits are far off-resonant from the main cav- 
ity mode, i.e. for Si, 62 3> 72. In this regime, both in the 
subradiant and in the non-subradiant scenarios, the main 
features of the dynamics can be obtained by looking at 
the effective dispersive Hamiltonian describing the cou- 
plin g o f the two qubits with a single-mode cavity field 
j22l. l23l. HH and remembering that this behavior must 
then be corrected taking into account the effect of the 
cavity losses. In Appendix B we derive the effective dis- 
persive Hamiltonian for this system, assuming that the 
cavity field is initially in the vacuum state, 

j—i 3 3 

(24) 

The first two terms in the Hamiltonian are propor- 
tional to crV ' cr_ and describe the Stark shifts due to 
the dispersive interaction with the cavity vacuum. The 
remaining terms describe an effective dipole-dipole cou- 
pling between the two atoms induced by the cavity mode. 
As we will see in the following these two terms play an 
essential role in the entanglement generation process. By 
looking at Eq. (f2"4"]) we notice that both the Stark shifts 
and the effective interaction strength between the qubits 
are now oc 72 2 /5i ] 2- 



In the dispersive regime the cavity is only virtually ex- 
cited, thus the photon loss is less important and the effec- 
tive decoherence rate due to the cavity decay is strongly 
suppressed to the advantage of the generation of entan- 
glement. As we will see in Sec. Ill A for the subradiant 
scenario, the effective decoherence rate due to the cavity 
decay in this case becomes (TZ 2 /S 2 )X. 



III. ENTANGLEMENT DYNAMICS 

To study the time evolution of the two-qubit entan- 
glement we use the concurrence C(t) [2l[. This is an 
entanglement measure related to the entanglement of for- 
mation, ranging from one for maximally entangled states 
to zero for separable ones. 

For the system of two qubits described by the reduced 
density matrix of Eq. {7J the concurrence takes a very 
simple form 

C(t) = 2\ Cl (t)c* 2 (t)\. (25) 

Such equation shows a relation between the behavior of 
the concurrence and the time evolution of the excitation 
shared by the two qubits. Having in mind the consider- 
ations of Sec. I D one may understand how, through a 
suitable choice of the detuning between the qubits and 
the cavity, it is possible to improve both the generation 
of entanglement and its preservation for long times. 

To better discuss the time evolution of the concurrence 
as a function of the initial amount of entanglement stored 
in the system, we consider a general initial states of the 
form given by Eq. (J5]) with 



coi = J—, c 02 = J ~y~ with -1<*<L 

Here, the separability parameter s is related to the initial 
concurrence as s 2 = 1 — C(0) 2 . 

Before describing in detail the dynamics in the sub- 
radiant scenario (Sec. Ill) and non-subradiant scenario 
(Sec. IV) , it is useful to recall the main features of the 
time evolution of the entanglement when u>\ = 0J2 = w c , 
i.e. in the resonant CcLSG, clS discussed in Ref. [H[. 

(i) The concurrence dynamics, as well as the value of 
the stationary concurrence, depends on the relative cou- 
pling strength, i.e. on the parameter r\. 

(ii) For certain entangled initial states, there exist at 
least one time instant t < 00 at which C(t) = 0, both in 
the strong coupling (good cavity) and weak coupling (bad 
cavity) limits, i.e. for A <C 72 and A 3> 72, respectively 
[See Fig. [2(a)]. 

(iii) In the weak coupling (bad cavity) limit, for an ini- 
tially factorized state, the reservoir creates entanglement 
and this is indicated by a monotonic increase in the value 
of the concurrence. 

(iv) In the weak coupling (bad cavity) limit, for an ini- 
tially entangled state, the reservoir causes entanglement 
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loss and the concurrence decreases with time until reach- 
ing, in some cases, the value zero, after which a small 
fraction of entanglement can be recreated [Sec Fig. [T] 
(a)]. 

(v) In the strong coupling (good cavity) limit, oscil- 
lations in the concurrence appear. For an initially fac- 
torized state there exist times at which the value of the 
concurrence is higher than the value of the stationary 
concurrence. 

In the next two sections we are going to study how the 
time evolution of the concurrence is modified in presence 
of detuning. 



IV. OFF-RESONANT ENTANGLEMENT IN 
THE SUBRADIANT SCENARIO 

We begin considering the case ui± = Whenever pos- 
sible, rather than discussing the exact expression of the 
concurrence, we will try to derive simpler approximated 
expressions which are useful for understanding the phys- 
ical processes taking place in the system. 



A. Bad cavity limit - Enhancement of the 
entanglement generation 

In the bad cavity case, e.g., for R = 1Z/X = 0.1, and 
for small values of the detuning S < 7Z, the behavior of 
the concurrence does not change appreciably compared 
to the resonant case. For values of the detuning 6 w TZ, 
i.e. when approaching the dispersive regime, the dynam- 
ics for an initially factorized state (s = 1) shows a mono- 
tonic increase towards the stationary value of the con- 
currence as in the resonant case. However, a significant 
change occurs in the bad cavity limit when the system 
is prepared in an initial entangled state. Indeed one can 
prove that in this regime, contrary to the resonant case, 
a finite time t such that C(t) = [See Fig. Q] (b)] does 
not exist anymore. 

We now focus on the dispersive regime i5 > A > K. 
If the qubit-system is initially entangled, e.g., for s = 0, 
the expression for the concurrence can be simplified as 
follows 

C(t) = \£\ w e~"?^" At , for n = 0, 1; (26) 
C(t) = \£\ 2 « e~ 2 ^ At , for n = (27) 

The equations above show that the concurrence vanishes 
with the decay rate (1Z 2 /S 2 )X when only one of the two 
qubits is coupled to the environment (ri = 0,1), and 
with 2(1Z 2 1 'S 2 )X when both qubits are identically coupled 
to the environment (for r\ = l/y/2). Since 1Z/5 <C 1 
this proves that in the dispersive regime the decay of 
entanglement is strongly inhibited compared to the res- 
onant regime since in this case the two atoms exchange 
energy only via the virtual excitation of the cavity field 



and therefore the cavity losses do not affect strongly the 
dynamics. 

For large enough detunings the entanglement shows 
oscillations as a function of time for all of the initial 
atomic states for which a finite stationary concurrence 
is obtained, C s ^ 0. Due to the presence of these oscilla- 
tions and for an initially factorized state, the concurrence 
reaches values greater than the stationary value C s even 
in the bad cavity limit, as shown in Fig. [2j For example, 
for n = n/3/2, R = 0.1 and 8 = 10A, at Xt « 2 x 10 3 
the concurrence reaches the value C = 0.92. For an ini- 
tially factorized state (s = 1) and for n = l/v2 we 
can derive the following approximated expression for the 
concurrence 

C(t) + e" 4 ^ At - 2e- 2 ^ At cos . (28) 

From this equation one sees that C(t) attains its max- 
imum value at t = imp- This formula also shows that 
the concurrence undergoes a series of damped oscillations 
with frequency 21Z 2 /S and decay rate 2{JZ/5) 2 X. 

With increasing detuning, the oscillations become 
more and more regular, quasi-periodic. The pattern is 
similar to the oscillations characterizing the strong cou- 
pling regime, but now the period is longer. As we will 
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FIG. 1: (Color online) Time evolution of the concurrence in 
the bad cavity limit (R = 0.1) with s = and = 0, for 
the cases of i) maximal stationary value, when r\ = \/3/2 
(black solid line), ii) symmetrical coupling n = l/y/2 (red 
dot-dashed line), and iii) only one coupled atom n = 0, 1 
(green dashed line). For each of such cases, we describe the 
entanglement dynamics in two different coupling regime: the 
resonant limit (left plot) and for 8% = 62 = 0.7A (right plot). 

see in Sec. Ill B, the generation of a high degree of en- 
tanglement in the dispersive regime for initially separable 
state can be achieved also in the good cavity limit. How- 
ever it is remarkable that already in the bad cavity limit, 
values of concurrence close to one can be generated. Our 
approach generalizes the results obtained for the disper- 
sive regime in Ref. [22| in the ideal cavity limit to the 
more realistic case of cavity losses. 

B. Good cavity limit - Entanglement quantum 
beats 

In the strong coupling case entanglement oscillations 
are present for any initial atomic state. Moreover, for 
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The corresponding cigcnencrgics arc given by 
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FIG. 2: (Color online) Time evolution of the concurrence in 
the bad cavity limit (R — 0.1) with s = 1, for the cases of 
i) maximal stationary value, n = v3/2 (black solid line), ii) 
symmetrical coupling n = l/v2 (red dot-dashed line), and 
iii) only one coupled atom ri = 0, 1 (green dashed line). All of 
the plots describe the dispersive regime with <5i = &2 = 10A. 



i" ~ A C K, when both atoms are effectively coupled to 
the cavity field, i.e., r% ^ 0,1, the dynamics of the con- 
currence is characterized by the occurrence of quantum 
beats, as shown in Fig. [3] For initially entangled states 
this phenomenon is more evident for <p = n because the 
value of stationary entanglement in this case is higher 
and the behavior of the concurrence is more regular. 

In order to better understand the origin of these entan- 
glement beats, we consider the case s = 1 and r\ = \j V2. 
For these values of the parameters, and for S ~ A <C 1Z, 
the expression of the concurrence can be written as fol- 
lows, 



C(t) 



1 



l + e 



-2Xt , 



cos(m) 4 - 2e- At cos(m) 2 cos(<5<). 

(29) 



The term 



cos(m) 2 cos(6t) = - cos(St) [1 + cos(2Kt)] 



in Eq. (|29|) . describing an oscillation at frequency 27?. 
modulated by a slower one with frequency 5, is responsi- 
ble for the occurrence of the quantum beats. 

To gain insight in the physical processes characterizing 
the dynamics, we consider the energy spectrum of the 
dressed states in the off-resonant case but in the absence 
of damping, as shown in Fig. 2] The diagonalization 
of the Tavis-Cummings Hamiltonian [See Eq. (j3"9")) in 
Appendix B] yields the dressed states 



I'M = 



1 



TZ 2 



+ w-\00)\l) R ), 



(30) 



1 



lu\ + 1Z 2 

\M = b/>->|o>fl. 



{-K\1>+)\0) r + uj + \00)\1) r ), 



(31) 
(32) 



U3± 

ujq = 5, 



= i (s ± \/m 2 



5 2 



(33) 
(34) 



where 1Z = ga^ is the vacuum Rabi frequency and 6 is 
the qubits-cavity detuning. 
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FIG. 3: (Color online) Time evolution of the concurrence in 
the good cavity limit (ii = 10) with s = 1, for the cases of 
i) maximal stationary value n = v3/2 (black solid line), ii) 
symmetrical coupling n = l/v2 (red dot-dashed line), and 
iii) only one coupled atom n = 0, 1 (green dashed line). The 
curves are drawn for small detuning, 6i = 82 = 0.7A; thus, 
outside the dispersive region. The inset shows the entangle- 
ment beat for the case i). 
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FIG. 4: (Color online) Energy spectrum of the dressed qubit- 
photon states in the case of small detuning (blue dashed line) 
and in the resonant coupling case (black solid line). 

On the other hand, the unperturbed states can be ex- 
pressed as a superposition of the |</>±),|</>o) eigenstates, 
with probability amplitudes evolving at frequencies io± 
and luq. The effect of the detuning is a shift of the qubits- 
cavity energy levels, thus the qubits-ficld coupling gives 
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FIG. 5: (Color online) Time evolution of the concurrence in 
the good cavity limit (R — 10), with s = and cj> = 0, for 
the cases of i) symmetrical coupling n = 1/V2 (black solid 
line), and ii) only one coupled atom n = 0, 1 (green dashed 
line). The two plots describe two different detuning regions: 
Si = 5 2 = 0.7A (left) and 8i= 5 2 = 50A (right). 



rise to a reversible energy exchange between unperturbed 
state at frequencies 272, 72 — 6/2 and 72. + 5/2. This is 
clearly seen, e.g., from the time evolution of the popula- 
tions 

|c 2 (t)| 2 = |(01|(0|e-^|01)|0) | 2 

r 4 / 5 

= r i + y I 1 + cos(272i)] + 2r\r\ cos (Tit) cos ( -t 

The equation above contains a term oscillating at fre- 
quency 27?., coming from the coupling between the 
dressed states |</>+) and | </>-), and a term oscillating at 



frequency 72 modulated by 8 coming from the interfer- 
ence between the oscillations at frequencies 1Z — 6/2 and 
1Z + 5/2 that couple the states \<fi+)-\4>o) and \4>-)-\4>q), 
respectively. 

In the discussion above we have disregarded the cavity 
losses. When they are taken into account one sees that 
the dressed energy splitting is resolved, and therefore the 
quantum beats will be visible, if 27?. is larger than the de- 
cay width A. This is achieved in strong coupling regime. 
Therefore, one does not observe quantum beats in the 
bad cavity case. 

We conclude this section studying how the detuning in- 
fluences the decay of entanglement, for an initially max- 
imally entangled state of the system, and the reservoir- 
induced entanglement generation, for an initial factorized 
state. 

When only one of the two qubits is effectively cou- 
pled to the cavity field, i.e. for r\ = 0, 1, for maximally 
entangled initial states (s = 0) in the resonant regime, 
5 = 0, the system performs damped oscillations between 
the states and \tp-}, which are equally populated at 
the beginning. Hence entanglement revivals with maxi- 
mum amplitude are present in the dynamics, as shown in 
Fig. [5] (a). Increasing the detuning, the amplitude of the 
oscillations decreases and the revivals disappear, while 
the frequency does not change appreciably, [See Fig. [5] 
(a)]. In this case the expression of the concurrence for 
small values of the detuning can be written as 



C{t) = \£\ « e- xt ' 2 \J cos{1Zt) 2 + P sin(7?i) 2 - A sin(72t) cos(72t), 

r 



(35) 



while for greater values of the detuning, the oscillations 
completely disappear and the concurrence decays expo- 
nentially 



C{t) = \£\ m 



-A/ 



(36) 



as shown in Fig. O (b). 

Finally, we note that, similarly to the behavior dis- 
cussed in the bad cavity limit, when the qubits arc ini- 
tially in a factorized state, the presence of the detuning 
enhances the generation of entanglement at short times 
compared to the resonant coupling illustrated in 

Fig. O In general, in the strongly dispersive regime, the 
qubits do not exchange energy with the cavity, which is 
only virtually excited. Thus a high degree of reservoir- 
induced entanglement can be generated both in the good 
and in the bad cavity limits. 



V. OFF-RESONANT ENTANGLEMENT IN THE 
NON-SUBRADIANT SCENARIO 

In this section, we analyze the more general situation 
in which the transition frequencies of the qubits are dif- 
ferent, u>i ^ 0J2, and both qubits are off-resonant with 
the cavity field. Due to the absence of a subradiant 
state, even a small value of the detunings 81,82 <C 1Z 
contributes to accelerate the decay of entanglement for 
every initial states. For an initially factorized state, in 
the bad cavity limit, the entanglement initially created 
via the interaction with the reservoir is rapidly destroyed 
as time evolves. In the good cavity limit entanglement 
oscillations are present and also quantum beats of entan- 
glement can be observed for 8\ , 62 ~ A <C 72.. 

We now consider in more detail the case in which the 
two qubits frequencies are symmetrically detuned from 
the central peak of the Lorentzian spectrum describing 
the field inside the cavity. In the dispersive region 5 3> 72, 
and for initially entangled states (s = 0), the concur- 
rence vanishes without manifesting a dominant depen- 
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FIG. 6: (Color online) Time evolution of the concurrence 
in the bad cavity limit (R = 0.1), with s = and cf> = 0, 
for the cases of i) maximal stationary value, corresponding 
to ri = y/3/2 (black solid line), ii) symmetrical coupling 
ri = l/y/2 (red dot-dashed line), and iii) only one coupled 
atom ri = 1 (green dashed line) and r\ = (blue dotted 
line). Two different detuning pairs are represented: the sym- 
metrical detuning with 8\ — — 0.7A, 82 = 0.7A (left plot) and 
the asymmetrical detuning with Si = — 0.5A, 82 = 0.9A (right 
plot). 



dence from r\ and <p. In other words, all the states 
initially entangled decay following the same behavior in 
such regime, as shown in Fig. [5] (a). This is in contrast 
to what we observed in all other regimes, where a de- 
pendence on the value of r\ is present. We stress once 
more that this feature seems to occur only for the case of 
symmetric detuning. Indeed, when introducing a small 
asymmetry in the value of the detunings the behavior of 
the concurrence shows again a dependence on the param- 
eter r\ , as illustrated in Fig. [5] (b) . 

In order to understand the peculiar behavior of the 
concurrence in the dispersive regime and for symmetric 
detunings we once more start by neglecting the effect of 
the cavity losses and use the dispersive Hamiltonian given 
by Eq. (pM|) . For symmetrical detunings §i = —62 this 
equation takes the form 

n K 2 rf (1) ( i)K 2 r 2 2 (2) (2) 

H eff = ^— o+gL' + — j- a+o-, (37) 

with S = \5 X \ = \S 2 \- 

Comparing Eq. (|37|) with Eq. (|24|) we notice that the 
terms describing the effective dipolc-dipolc coupling in- 
duced by the cavity mode are here absent. Therefore the 
only remaining effect is the entanglement decay induced 
by the cavity losses. The decay rate, however, does not 
depend on the relative coupling parameter t\ but only 
on the total coupling strength olt via the vacuum Rabi 
frequency 1Z. This explains why, even when the cavity 
losses arc taken into account, the time evolution of the 
concurrence for symmetric detunings does not depend on 
r\. When a small asymmetry in the detunings is intro- 
duced, the dipole-dipole effective coupling terms are non- 
zero and, due to the presence of T\ and r-z in the effective 
dipole-dipole coupling strength, the dynamics becomes 
again dependent on r\ . 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have provided a complete analysis of 
the exact dynamics of the entanglement for two qubits 
interacting with a common zero-temperature reservoir in 
the off-resonant case. We have presented a general ana- 
lytical solution for the two-qubit dynamics without per- 
forming the Born-Markov approximation. In the case of 
a Lorentzian spectrum, describing, e.g., the electromag- 
netic field inside a single mode lossy cavity, we have ob- 
tained explicit expressions for the reduced density matrix 
and for the concurrence. The availability of the exact so- 
lution allowed us to look at the entanglement dynamics 
both in the weak coupling (bad cavity) and in the strong 
coupling (good cavity) limits. 

If the two qubits are initially disentangled, the interac- 
tion with the common reservoir generates entanglement. 
Our results demonstrate that a high degree of entangle- 
ment can be generated in this way, especially in the dis- 
persive regime, and even in the bad cavity limit. For ini- 
tially entangled states, the concurrence decay is slowed 
down when the qubits are detuned from the peak of the 
Lorentzian. In this case, indeed, the cavity losses affect 
less the atoms dynamics since the effective atom-atom 
interaction is mediated by virtual photon exchange. 

In general, the entanglement dynamics is strongly sen- 
sitive to the relative coupling parameter n, indicating 
how strongly each of the two qubit is individually coupled 
to the e.m. field. Only when the qubits frequencies are 
symmetrically detuned from the main cavity frequency, 
in the dispersive regime, the dependence on the relative 
coupling disappears. Finally we have discovered that, 
in the strong coupling regime, for intermediate values of 
the detuning, the dynamics of the concurrence shows the 
occurrence of quantum beats. We have given a phys- 
ical interpretation of this phenomenon in terms of the 
quantum interference between the transitions among the 
dressed states of the atomic system. 

We believe that our results contribute in shedding light 
on the behavior of quantum entanglement in realistic con- 
ditions, that is when the effect of the environment on the 
quantum system is taken into account. For this reason 
they have both a fundamental and an applicative value 
and they indicate how rich the dynamics of this system 
can be. The model we have studied can be employed to 
describe both trapped ions in optical cavities [251 ] and cir- 
cuit cavity QED dynamics [26|, [U [28| ■ In both physical 
contexts, the observation of the effects we have discussed 
should be achievable with the current experimental tech- 
nologies. 
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Appendix A - Analytic solution for the probability 
amplitudes 

In this Appendix we briefly discuss the structure of 
the analytical solutions of Eqs. (fll3|) - (fl"4"]) for the proba- 
bility amplitudes ci j2 (i) and how they can be obtained 
applying the Laplace transform method. We note that 
the solutions obtained in this way are exact since we do 
not perform any kind of approximation. 

The solution of the Laplace transformed amplitudes 
ci, 2 (s), obtained from Eqs. (fl~5|) -(fT6 |) can be written as 
the sum of three ratios having denominators (s — s;), 
where Si are the roots of the cubic equation. 

s 3 + A r s 2 + B jS + Cj = 0, (j = 1, 2) (38) 

where 



can apply the canonical transformation defined by the 
unitary operator [23|, [24| 

e aS = e »i y + - ) (39) 

with IZrj = goij . This procedure is correct to the second 
order in the coupling to the cavity, and, limiting our- 
selves to this approximation, we can write the effective 
Hamiltonian as follows 

2 

H eff = e aS He~ aS ~H + a[S,H} + ^-[S, [S, H}]. 

Assuming that the cavity field is initially in the vacuum 
state, H e ff takes the form 

2 n 2 r 2 K 2 nr 2 / (1) (2) (2) 

i=l J J 

(40) 

(i) (i) 

where the terms proportional to aY a_ describe the 
Stark shifts due to the dispersive interaction, while the 
last two terms describe the dipole-dipole coupling be- 
tween the two atoms induced by the cavity mode through 
the exchange of virtual cavity photons. 



Aifl = X + i (5 2 ,-i -2<5i. 2 ), 

Bifl = II 2 - Sf t2 + 5 l 5 2 +i fax - Si, 2 ) A, 

d,2 = %H 2 r 2 2 (5 2; i-5i, 2 ). 

The amplitudes Ci t2 (t), obtained by inverse Laplace 
transform will then be the sum of three damped oscil- 
lating terms having, in general, a complicated structure. 
Only for the case lo\ = lo 2 a simple analytical expressions 
for the probability amplitudes can be obtained, whereas 
in the general case there is no simple solution. This is 
because, when lo\ = cu 2 the cubic equation can be writ- 
ten as a product of polynomials of first and second order 
having always one root coincident with zero. In this case 
one can write the amplitudes in the simple form given by 
Eqs. CII-CEII). 



Appendix B - Effective dispersive Hamiltonian 

The Hamiltonian describing the interaction between 
two-qubit systems and the quantized cavity mode is given 
by 

2 

+ \g (aia^ + a 2 a { + ] ^ b + h.c. . 

To obtain the effective Hamiltonian describing the in- 
teraction with the cavity in the dispersive regime, one 



Appendix C - Approximate expressions of the 
concurrence 

In this Appendix we derive approximate expressions 
for the amplitudes ci, 2 (t) in the case of large (and equal) 
detuning 5 ^> 1Z. 

For this purpose, we expand the term Q = 
^/A 2 — Sl^ — i25\ as follows, 

n .,( 1 _^). < (, + 2|!). m 

The temporal evolution described by £ (t) can then be 
written as 







cosh | — ) 




v 2; 





„ e -f?(A+«)*. 

For the sake of simplicity we consider here the case 
s = 1 and n = 1/ \/2. However the time evolution of the 
concurrence has features in common with all of the other 
cases: 

C(t) = 2|c i(f)||c 2 (t)| 

= y(l + \E(t)\ 2 f-(2Re[E(t)]f 

-l^ 1+ .-V» (4,). 
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On other hand, for small detunings of the order of A, Therefore, for the case s = 1 and n = l/v^ the time 
outside the dispersive region i « R, the approximate evolution of the concurrence is given by 
form of f2 is given by 



A S 

il w — - i2K, 
2TI 



(42) 



so that the time evolution is described by the function 



£{t) 



-(\-iS)t/2 



-(X-iS)t/2 



cosh I — - ) — 



S + iX . 



sinh — — 



Of 



c(t) = -#TwFWiF 



-■sjl + e- 2Xt cos(m) 4 - 2e- xt cos(Kt) 2 cos(St). 



2K 

X s 

cos (Tit) sin (Tit) +i sin (Tit) 

y ' 2TI y ' 2TI 
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